Hyperfine-induced valley mixing and the spin-valley blockade in carbon-based 

quantum dots 
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Hyperfine interaction (HFI) in carbon nanotube and graphene quantum dots is due to the presence 
of ^''C atoms. We theoretically show that in these structures the short-range nature of the HFI 
gives rise to a coupling between the valley degree of freedom of the electron and the nuclear spin, 
in addition to the usual electron spin-nuclear spin coupling. We predict that this property of the 
HFI affects the Pauli blockade transport in carbon-based double quantum dots. In particular, we 
show that transport is blocked only if both the spin and the valley degeneracies of the quantum dot 
levels are lifted, e.g., by an appropriately oriented magnetic field. The blockade is caused by four 
"supertriplet" states in the (1,1) charge configuration. 

PACS numbers: 73.63.Kv, 73.63.Fg, 73.23.Hk, 31.30.Gs 



In the past decade, fundamental steps have been 
made towards the realization of quantum information 
processing, including isolation, manipulation, and read- 
out of single electron spins in the solid state.— How- 
ever, the majority of the existing quantum dot (QD) 
spin qubits is fabricated in material systems, where 
hyperfine interaction (HFI) with nuclear spins limits 
the device performance via spin decoherence. Car- 
bon structures, such as carbon nanotubes (CNTs) or 
graphene, are expected to have weak HFI, due to the 
small 1% natural abundance of spin-carrying -'^'^C nu- 
clei. This expectation has motivated intensive theoret- 
ical investigation2ii^i^>^ and the experimental realization 
of QDs in carbon nanostructures-^i^i^iiaiiiii^ii^iiii^ A 
further perspective of carbon-based quantum information 
processing has been opened by proposals suggesting to 
utilize the valley degree of freedom of the delocalized elec- 
trons as a qubit ji^iii Relaxation and decoherence mech- 
anisms of these valley-qubits are yet to be explored. One 
possible source of those is short-range disorder, which is 
known to couple the two different valley states.— 

Double quantum dot (DQD) structures in the two- 
electron regime are particularly well suited for study- 
ing the effects of HFI .i^i^o 121^.23 ^-j^^ so-called Pauli 
blockade regime, the measurement of the direct current 
through a serially coupled GaAs DQD, as a function of 
the external magnetic field, has been used to infer the 
hyperfine energy scale . ^"'^"'^ Another experiment in GaAs 
DQDs showed that HFI can be utilized to perform coher- 
ent rotations between the states of a singlet-triplet qubit 
and, at the same time, acts as a source of decoherence.— 

The effect of the HFI in carbon is most pronounced 
in fully ^■^C-enriched samples. ^•^•^^'^'^'-^'^'^^ Such nanotube 
DQD devices have been used recently to estimate the 
energy scale of the atomic HFI as ~ 100/ieV, using 
transport^'' and singlet-triplet dephasing timc^ measure- 
ments. In contrast, theory predicts an atomic hyperfine 
energy scale ~ l^eVi^i^i This discrepancy between the- 
ory and experiment, together with unexplained features 
of the dephasing time measurements of RefH, show that 
additional theoretical efforts have to be made to gain 



a complete understanding of the role of HFI in carbon- 
based QDs. 

Here, we study the influence of the ^'^C nuclear spins 
on the spin and valley degrees of freedom of the electrons 
in carbon-based QDs. In particular, we derive the 4 x 
4 Hamiltonian describing the effect of HFI on a single 
fourfold (spin and valley) degenerate QD energy level. 
We find that due to the short-range nature of the HFI, it 
couples the nuclear spins not only with the spin, but also 
with the valley degree of freedom of the electron. The 
effective hyperfine Hamiltonian can be expressed as 
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Here S = ( Sx,Sy,Sz)/2 is the spin operator, tq is the 
unit operator in valley space, Si (r^) denotes the Pauli 
matrices acting in spin (valley) space, and the quanti- 
ties are different linear combinations of the in- 
dividual nuclear spin operators (see below) . Equation ([1]) 
should be contrasted with the widely used 2x2 hyper- 
fine Hamiltonian i/hf.GaAs = S ■ h, which describes the 
effect of the nuclear spin on a twofold degenerate level in 
a GaAs QD, and incorporates only a single Overhauser 
field h. We estimate that the order of magnitude of the 
valley-conserving (^ tq) and valley-mixing parts of TJhf 
are the same. 

As a physical consequence of the valley coupling due 
to the HFI, we predict that the response of the Pauli 
blockade leakage current through a carbon-based DQD 
to an applied external magnetic field is remarkably dif- 
ferent from the case of GaAs DQDs. In the Pauli or 
spin blockade regime^i transport from the source to the 
drain through a serially coupled DQD occurs via the 
(0,1) ^ (1,1) ^ (0,2) ^ (0,1) cycle, {ul^ur) denot- 
ing the charge state with [nR] electrons in the left 
(right) QD [see Fig. □ (inset)]. In a GaAs DQD, block- 
ing of the current occurs when there is at least one two- 
electron energy eigenstate in the (1,1) charge configu- 
ration having a spin wave function, which is symmetric 
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under particle exchange (i.e., a triplet). Due to HFI, this 
condition is achieved only in the presence of an external 
magnetic field, which splits two triplet states apart from 
the singlet state and prevents hyperfine-induced mixing 
of thosei^iSi In carbon-based QDs with fourfold level de- 
generacy, the current is blocked only if there is at least 
one energy eigenstate in the (1,1) charge configuration 
with a combined spin-valley wave function which is sym- 
metric under particle exchange ('supertriplet'). In order 
to distinguish this effect from the spin blockade in con- 
ventional DQDs, we call it the spin-valley blockade. We 
show that in contrast to GaAs, in carbon DQDs a spin 
(Zeeman) splitting is insufficient to maintain the block- 
ade in the presence of HFI, which, however, can be re- 
covered by simultaneously introducing spin and valley 
splittings. 

We consider the lowest-lying orbital level of an elec- 
trostatically defined QD in monolayer graphene with a 
finite gap.'^ (Nevertheless, the concepts and the formal- 
ism we use are readily generalizable to other types of 
graphene or CNT QDs.) In the absence of nuclear spins 
and external magnetic field, this level is fourfold (spin 
and valley) degenerate. In the presence of an external 
magnetic field, both the spin and valley degeneracy can 
be split: an in-plane magnetic field causes only a spin 
splitting As = QefJ-sB via the Zeeman effect, whereas 
an out-of-plane (z) component introduces a valley split- 
ting Ay{Bz) as welli^ We consider the regime, where 
the energy difference between the lowest-lying and the 
second orbital levels is much larger than As and A^,. 
We describe the system using the tight-binding (TB) 
model. The TB wave functions, corresponding to the 
four sublevels of the lowest-lying orbital level and char- 
acterized by the spin and valley quantum numbers s G 
(t,i) = (+,-) and V e {K,K') = (+,-), are {i^sv)ia = 
V^ie'''^"'''-^a\ria)Xs- Here a G {A,B} is the sub- 
lattice index, I is the unit cell index, Slccii is the unit cell 
area, ri^r is the position of the carbon atom on sublattice 
a in the Ith unit cell, and x+ — (1, 0) and x- = (0; 1) are 
the two possible spin states. The four smoothly varying 
functions '^a can be obtained by solving the Dirac-like 
envelope function equation. '^'^^ The functions 5'^"^ and 
are normalized: /dV (|*^''^(r)|2 + |*^''^(r)|2) = 1 

and 'Zla \i'^^sv)la\^ = 1- 

The nuclear spin of the carbon atom on site la is de- 
noted by I la, being zero if the atom is a ■^'^C and a 
spin- 1/2 operator if the atom is ^^C. In graphene and 
CNTs, the HFI is known to be short-range: its major 
contribution to the TB Hamiltonian is the on-site ma- 
trix element on the site of the nuclear spini^ Therefore 
the elements of the TB Hamiltonian matrix, describ- 
ing the HFI, are (i?hf,tb);<T,iv = SwS^^' SAIi^, where 
A = dia,g{Ax,Ay, A^) is a diagonal matrix. 

The effective Hamiltonian describing the influence of 
the HFI on the four sublevels is constructed by express- 
ing -ffhf.tb in the subspace spanned by the four TB wave 
functions ipsv The resulting matrix can be expressed in 




FIG. 1; (Color online) Averaged leakage current (7) due to 
hyperfine interaction in a graphene DQD, as a function of 
spin and valley splittings, /max ~ 0.34erin. Inset: transport 
through a DQD. Dot R is always occupied by at least one 
electron. 



terms of the spin and valley operators as shown in Eq. 
(P) using 
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^'CT^''*(f ;o-)^'i \ria)- Note that in the presence of time 
reversal symmetry /i*-^-* = 0. Under normal conditions, 
the nuclear spins are completely randomized by thermal 
fluctuations, which implies that the components of the 
hyperfine fields have zero mean and their variances are 
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We have used {{ha)"]) = v/A and v denotes the 
abundance of -'^'^C atoms in the QD. Using the slowly 
varying nature of the envelope functions, we find 

Ez. (^Z;'"')' = lUia \9ia?- In the TV - ex. limit [N 
is the number of atoms in the QD) the distributions of 
the hyperfine fields /i^'^'' converge to Gaussians. Further- 



more, h^^^ and /i^-^"* become independent from hf\ 



and each other. 

To determine the variances of the hyperfine field pre- 

(v) 

cisely, one needs to know the envelope functions '^a ■ 
However, a simple estimation can be given using the as- 
sumption ~ \/^lcc\\N. This choice satisfies the 

normalization condition, and it results in {{h^P^ ) = 
= A]v/m and {(hf ^)^) = 0. We 



A]i,/AN, ((/if 

will use these values in the following calculations. 

The Hamiltonian H\^f in Eq. ([1]) shows that the HFI in 
carbon-based QDs results both in spin and valley mix- 
ing, and our estimates of the hyperfine field variances 



3 



suggest that the valley-conserving and valley- mixing con- 
tributions in 7?hf have the same order of magnitude. We 
emphasize that this result is due to the short-range na- 
ture of the HFI. We anticipate a similar result in sili- 
con QDs,^'*'^^ where electrons also possess a valley degree 
of freedom. We highlight three possible physical conse- 
quences of the valley- mixing nature of -ffhf- (i) HFI can 
be a source of valley relaxation and decoherence in 'val- 
leytronics' devices, such as valley filters, valley valves^ or 
QD valley qubitsJ^ (ii) Spin-orbit couphng in CNT QDs 
can split the fourfold degenerate ground state dot level 
into two doublets {K^, K' [ and if't, K [) at B = 0,1 
Spin-orbit coupling is stronger than the HFI; hence, a 
fully valley-conserving iJhf would not be able to cause 
mixing within or between the Kramers doublets. Accord- 
ing to our result, _ffhf is not valley-conserving and causes 
mixing within the Kramers doublet. This finding may 
facilitate the understanding of yet unexplained features 
of recent experiments on ^'^C nanotubes^iiP. (iii) As we 
show below, the valley-mixing character of the HFI in- 
troduces remarkable features in the behavior of the Pauli 
blockade leakage current through a carbon-based DQD. 

Our goal is to calculate the average leakage current 
through a DQD in the Pauli blockade regime as a function 
of spin and valley splittings. To this end, we have gen- 
eralized the master equation formalism used in Ref 21. 
We focus on the parameter regime where the result Eq. 
(11) of Ref Hi] is valid: (i) The energy detuning A = 
£;(0, 2) -£;(!, 1) between the (0,2) and (1,1) charge states 
is much larger than the coherent tunneling strength be- 
tween the two dots (t) , the characteristic energy scale de- 
scribing the HFI [/ihf = ^{v/AN) J2j Af\ , and the spin 

and valley splittings (Ag, A„). This condition suppresses 
the coherent tunneling between the dots; therefore, the 
hybridization between (1,1) and (0,2) charge states is 
negligible, (ii) The (1,1) (0,2) transition is an ener- 
getically downhill inelastic tunneling process, character- 
ized by the rate Pin. (iii) i^/|A| <C /ihf, which enables 
one to neglect the exchange splitting within (1,1) charge 
states, (iv) Pl » P^j ^ Pin, where the rate Pl [Pi?] 
describes the (0,1) (1,1) [(0,2) (0,1)] transition, 
(v) Pi < ft,hf , which enables one to use a classical master 
equation to describe the transport process. Studying the 
Pauli blockade problem under the above specified condi- 
tions is motivated by the fact that the theory for GaAs 
DQDs in this parameter regime has been successful in de- 
scribing recent experimental results:™ We also adopt the 
constant interaction feature of the model used in Ref 21. 
Investigation of parameter regimes where the Coulomb 
interaction results in Wigner molecule formation^2i^i22 
is beyond the scope of this Rapid Communication. 

Our approach deviates from the one used in Refl2ll 
in the form of the single-particle Hamiltonians of the 
two QDs, incorporating external magnetic field and 
HFI: Hd = -/^D.magn + Hd,m, where D e {L,R}, 
Ho.magn = gc^^BSo ' B + A„(B^)t2,d/2 and i?D,hf = 

■ Yji=Q.x,y,z^^D'^i,D- According to our above analy- 
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FIG. 2: (Color online) Averaged leakage current (/) as a func- 
tion of (upper solid line, black) spin splitting As for zero valley 
splitting A„ = 0, (lower solid line, green) valley splitting A„ 
for As = 9, (dashed) As in a GaAs DQD [dashed line based 
on Eq. (11) of RefiJ. 



sis, we treat hj^' as uncorrelated stationary random 
fields^i and use our previous estimates for their variances. 
For our purposes, we can neglect the spin-anisotropy of 
the HFI,^'^'* and use Aj = Aiso and therefore /ihf = 
A\so\/'iv / AN . Evaluating /ihf with the measured^iiiS 
value ^iso = 100//eV and assuming iV = 7.5 X 10^ atoms 
in a QD we find /ihf ~ 30neV (SOOneV) for abun- 
dance V — 1% (100%); using the theoretical estimate^'— 
Also = l/xeV, we find /ihf ~ 0.3neV (3neV). Under the 
above specified conditions, we set up a classical mas- 
ter equation for the occupation probabilities of the (1,1) 
eigenstates of Hl + Hr. By solving the master equa- 
tion numerically, we calculate the stationary current / 
for many realizations of the nuclear fields and average 
those to obtain (/). In Fig. [1] we plot (/) as a function 
of spin and valley splittings obtained by averaging for 
200 random realizations of the hyperfine fields. Two cuts 
along the horizontal (black) and vertical (green) lines of 
Fig. [1] are shown on Fig. [21 obtained by averaging over 
3000 random realizations of the nuclear fields. For com- 
parison, in Fig. [2]we also plotted the result corresponding 
to the case of GaAs DQDs (dashed). 

One of the characteristic features of the spin-valley 
blockade is the cross-shaped pattern in the (/(As,A^)) 
density plot in Fig. [TJ This pattern indicates that the 
current is strongly suppressed only if both the spin and 
valley splittings (A^ and A^,) exceed the energy scale of 
the hyperfine coupling. As mentioned earlier, a physi- 
cal situation where As is finite but A^ = 0, is when an 
in-plane magnetic field is applied to the system. This 
situation corresponds to the black line in Fig. [T]and the 
solid black curve (/(A^)) in Fig. [2l It is apparent from 
Fig. [2] that the current decreases with increasing A^; but 
instead of approaching zero, it saturates to a finite value 
around 0.2ePin. This is in stark contrast to the case of 
GaAs DQDs, where As > /ihf leads to a sharp decay of 
the current (dashed curve). 

To present a qualitative interpretation of our results. 
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FIG. 3: (Color online) Schematic energy diagrams of the (1,1) 
charge configuration corresponding to the three highlighted 
points of Fig. [1] The effect of hyperfine interaction is ex- 
cluded. Dashed lines; transport-blocking supertriplet states. 



we make use of the analogy between the spin and valley 
degrees of freedom. We describe the states of the (1,1) 
charge configuration by using the simultaneous eigenbasis 
of the total spin operator {Sl + Sr)^, the spin projection 
on the direction of the magnetic field (Sl+Sr)-B/B, the 
total valley operator (x^ + t^)^/4, and the z-componcnt 
of the valley operator (r^x + Tz,r)/2. The correspond- 
ing quantum numbers are s £ {0, 1}, nis S {—s, . . . , s}, 
V € {0, 1}, and € {—v^ . . . We denote these ba- 
sis states with \s,ms,v^mv) ■ These are eigenstates of 
the system Hamiltonian in the absence of HFL The com- 
bined spin-valley wave functions of the ten states fulfill- 
ing s = V are supertriplets; therefore, these states cannot 
be squeezed into a (0,2) charge configuration. In con- 
trast, the spin- valley wave functions of the six states with 
s ^ V are supersinglets; hence, their transition to (0,2) is 
allowed. 

The energy diagram of the 16 states of the (1,1) 
charge configuration, corresponding to the three high- 
lighted points of Fig. [1] are presented in Fig. [31 Figure 
[5^ shows the situation, where = = 0. In this case, 
there is a 16-fold degenerate level, and the HFI mixes 
supersinglet and supertriplet states effectively. This re- 
sults in a maximal current through the DQD. The Zee- 
man effect splits the states with different quantum 



numbers (Fig. [Sja) and suppresses hyperfine-induced hy- 
bridization between them if A<, > /ihf, which leads to 
a decrease in the leakage current. However, the valley 
mixing contribution of the HFI still induces strong mix- 
ing within the states with the same nis- This mixing 
prevents the appearance of "pure" supertriplet energy 
eigenstates which would block the transport; therefore, 
the current does not drop to zero. As mentioned earlier, 
this behavior is in contrast to the case of GaAs DQDs. 
Figure ^ shows the energy diagram when both Ag and 
A„ are finite. If those are larger than the HFI, then 
the four supertriplet states |1,±1,1,±1) become decou- 
pled from supersinglets. Thus, the system gets trapped 
whenever any of these four states is occupied during the 
transport process, which results in a strong suppression 
of the current. Note that only two blocked states remain 
if As w A„. In that case, |1, 1, 1, -1) and |1,-1, 1, 1) be- 
come degenerate with the fourfold degenerate |s,0, w,0) 
and mix with those due to HFI, slightly enhancing the 
current, visible along the diagonal |As| = |A„| lines in 

Fig.m 

Another characteristic of the spin-valley blockade is 
the appearance of a dip in the green (/(A„)) curve in 
Fig. [2] at A^ = 0. Similar dip structures have been 
predicted^i^i and measure d^°i^°'^^ in conventional semi- 
conductors and they were attributed to various micro- 
scopic origins, including cotunneling, spin-orbit interac- 
tion, and exchange coupling. In our case, the dip has a 
different origin: it is due to the strong valley anisotropy 
of the HFI, i.e., that in our above estimations van- 
ishes and therefore i/hf does not include the operator. 

We have established the form of the Hamiltonian de- 
scribing the effect of HFI on a fourfold degenerate energy 
level in a carbon-based QD. We have found that the short 
range nature of the HFI leads to a significant nuclear 
spin-electron valley coupling. We have calculated the ef- 
fect of this interaction on the leakage current through 
a DQD in the Pauli blockade regime. Our findings may 
have profound consequences for both spin and valley ma- 
nipulation in carbon-based QDs. 
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